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ABSTRACT 

FU Tau belongs to a rare class of young, wide brown dwarf binaries. We have resolved the 
system in a Chandra X-ray observation and detected only the primary, FU Tau A. Hard X-ray 
emission, presumably from a corona, is present but, unexpectedly, we detect also a strong and 
unusually soft component from FU Tau A. Its X-ray properties, so far unique among brown 
dwarfs, are very similar to those of the T Tauri star TW Hya. The analogy with TW Hya sug- 
gests that the dominating soft X-ray component can be explained by emission from accretion 
shocks. However, the typical free-fall velocities of a brown dwarf are too low for an interpre- 
tation of the observed X-ray temperature as post-shock region. On the other hand, velocities 
in excess of the free-fall speed are derived from archival optical spectroscopy, and indepen- 
dent pieces of evidence for strong accretion in FU Tau A are found in optical photometry. 
The high X-ray luminosity of FU Tau A coincides with a high bolometric luminosity con- 
firming an unexplained trend among young brown dwarfs. In fact, FUTau A is overluminous 
with respect to evolutionary models while FU Tau B is on the 1 Myr isochrone suggesting 
non-contemporaneous formation of the two components in the binary. The extreme youth of 
FU Tau A could be responsible for its peculiar X-ray properties, in terms of atypical magnetic 
activity or accretion. Alternatively, rotation and magnetic field effects may reduce the effi- 
ciency of convection which in turn affects the effective temperature and radius of FU Tau A 
shifting its position in the HR diagram. Although there is no direct prove of this latter scenario 
so far we present arguments for its plausibility. 
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1 INTRODUCTION sured from infrared (IR) photometry to be ~ 2 mag for the primary 

and < 1 mag for the secondary. The spectral energy distributions 
The youna binary brown dwarf (BD) FUTau has aroused inter- r> . i ■ ■ i- i r ■ n 

ot both binary components indicate the presence or circumstellar 
est in the star formation community for two particular aspects j- i i ■ .i i , i • j i ■ n ■ 

^ ^ ^ disks by excess emission over the photospheric model m the Spitzer 

both having important implications for BD formation theories (see j m i_ j t jj-.- • i i r T-xTrr. a . 

I — ; — " II r I *^ r ■ mid-lR bands. In addition, there IS a blue excess tor FU Tau A that 

ILuhman et aIJI009af ): (1) It is one of only a handful of known bi- 



nary BDs with wide separation (800 AU). (2) It is located in the 



might indicate ongoing accretion. Both FU Tau A and B have strong 

Ha emission in low-resolution spectra (equivalent widths of 93 and 
B215 dark cloud in Taurus that hosts only one other T Tauri star _„ s i ^ ■ i-i i ■ , 

70 A, respectively), again a likely accretion signature. 

(TTS). It IS the only young (< 5 Myr) BD known in such a re- 
mote location and, therefore, it has formed most likely in isolation Widely separated BDs are important calibrators of sub-stellar 
and not in a stellar c luster or aggregate. From low-resolution opti- magnetic activity because they allow to examine the X-ray char- 
cal spectra lUuhman e t al. (009a) measured spectral types of M7.25 acteristics of two presumably coeval targets with different (sub- 
and M9.25 for the two components FUTau A and FUTauB. They )stellar properties (T^ff, Lboi, etc.) in a critical parameter range, 
derive masses of 0.05 and 0.015 Mq, respectively, comparing Magnetic activity (as probed by Ha and X- ray luminosity) seems 
the position of the two BDs in the HR dia gram to pre-MS evolu- to decay strongly at late-M spectral types iMohantv et alj|200j : 
tionary models bv' Baraffe et al.l (1 199 8) and C habrier et all bOOd) . IStelzer et al]|200d ). Chandra is the only X-ray instrument capable 
They find that FU Tau A is significantly overluminous relative to the of resolving FU Tau and similar BD binaries. In an earlier study we 
youngest (1 Myr) isochrone of the models. The extinction was mea- have resolved the wide BD binary in the C hamaeleon star-formin g 

region, 2MASS II 101 1926-7732383 AB ( Stelzer & Micelall2007h . 

In a continuation of our search for coronal activity in sub-stellar 
* E-mail: stelzer@astropa.inaf.it twins we have obtained a Chandra X-ray observation of FUTau. 
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Figure 1. Chandra/ ACIS image of the FUTau binary. 2 MASS positions 
are marked with crosses. The secondary is not detected. 



Here we combine the analysis of the X-ray data with archival op- 
tical spectroscopy and photometry. The primary turns out to have 
unexpected X-ray properties. This observation shows how incom- 
plete our knowledge of the high-energy emission of BDs is. 



2 DATA ANALYSIS AND RESULTS 
2.1 X-ray observation 

On 22 Oct 2009 FUTau was observed for 55ksec with 
Chandra/AClS-S (Obs-ID 10984). The data analysis was carried 
out with the CIAO packagfl version 4. 1 , we started with the level 1 
events file provided by the Chandra X-ray Center. The steps per- 
forme d to convert these data t o an event level 2 file are described 
e.g. in lStelzer & Micelal ( l2007h . 

The X-ray image is almost devoid of X-ray sources, not un- 
expectedly as high extinction through the B 215 cloud obscures the 
background sky. However, bright X-ray emission is found to co- 
incide with the 2 MASS position of FUTau A. For detecting the 
X-ray source(s) associated with the two components of the BD bi- 
nary, source detection was restricted to a 100 x 100 pixels wide 
image (1 pixel — 0.5") and a congruent, monochromatic exposure 
map for 1.5 keV centered on the 2 MASS position of FUTau A. 
Source detection was carried out with the WAVDETECT algorithm 
with wavelet scales between 1 and 8 in steps of \/2- We tested a 
range of detection significance thresholds, and found that ath = 
10~^ avoided spurious detections and at the same time separated 
close emission components. The X-ray image in the region around 
FUTau A is shown in Fig.[T] There is only one X-ray source in the 
surroundings of the BD binary. It is associated with FUTau A. The 
companion FU Tau B remains undetected. 

We calculated the source count rates in the following way: A 

^ CIAO is made available by the Chandra X-ray Center and can be down- 
loaded from |http://cxc.harvard.edu/ciao/download/| 



Table 1. X-ray parameters of the FU Tau binary. 

Object Off ax Counts Expo log Lx \og{-^^-) 
['] in0.3-8keV [sec] [erg/s] 

FUTau A 029 603.8 ± 24.6 54153 207 -3.2 

FUTauB 0.39 < 4.7 54053 < 27.2 < -3.8 
Bolometric luminosities from^uhmanetaf ^009a). 

circular source photon extraction region was defined as the area 
that contains 95 % of the point spread function at the position 
of FU Tau A. The background was extracted individually from a 
squared region centered on the source extraction area and several 
times larger than the latter one. A circular area centered on the po- 
sition of the X-ray source was excluded from the background area. 
The S/N was computed from the counts summed in the source and 
background areas, respectively, after applying the appropriate area 
scaling factor to the background counts. In practice, the background 
is very low (a fraction of a count in the source extraction area). Fi- 
nally, the count rate of FU Tau A was obtained using the exposure 
time at the source position extracted from the exposure map. We 
have estimated a 95 % confidence upper limit for the coun t rate a t 
the position of FU Tau B using the algorithm of I Kraft et alj ( Il99 ih . 

In Table[T]we summarize some X-ray parameters of the binary. 
The distance assumed for the estimate of luminosities is 140 pc. For 
FU Tau A we report the value of Lx extracted from its spectrum (see 
below), while for FU Tau B we used PIMM^I to convert the count 
rate upper limit to a flux limit. This conversion depends on the un- 
known spectral parameters. The value given in Table [prefers to an 
assumed A'^h ~ 2 ■ 10^^ cni~^ (corresponding to the upper limit 
measured for the A-v of FUTauB) and kT — 1 keV (as typical for 
a coronal plasma). Note that the influence of TVh on the count-to- 
flux conversion is much stronger than that of the temperature. 

For FU Tau A a lightcurve was extracted and searched for vari- 
ability with a maximum likelihood method that di vides the se- 
quenc e of photons in intervals of constant signal (see IStelzer et alJ 
|2007|) and, independently, with the Kolmogorov-Smimov (KS) 
test. There is no significant variability in the X-ray count rate of 
FUTau A (KS-test probability is 0.34), such that the bright emis- 
sion can not be attributed to a flare. 

The high signal has allowed us to fit the X-ray spectrum of 
FUTau A. An individual response matrix and auxiliary response 
were extracted for the position of FU Tau A using standard CIAO 
tools. As mentioned above, the background of ACIS is negligibly 
low. We fitted the spectrum, rebinned to a minimum of 15 counts 
per bin, in the XSPEC 12.5.0 environment with a one- or two- 
temperature thermal model subject to absorption. The spectrum of 
FU Tau A can be described by an absorbed 2-T thermal model (see 
Tabl e [2] for the best fit param eters). For standard extinction laws 
(e.g. lPredehl & SchmittI 19951) the observed X-ray absorption, A^h, 
is compatible with the extinction of Ay — 2 mag measured in the 
IR. 

In summary, with more than 600 counts the Chandra obser- 
vation of FU Tau A represents the highest quality X-ray data ob- 
tained so far for a (young) BD. The most remarkable fact is that the 
emission is dominated by a cool plasma component of temperature 
kT — 0.24 keV. The soft plasma has four times more emission 
measure (EM) than the hotter component that has a typical coronal 

^ The Portable Interactive Multi-Mission Simulator (PIMMS) is accessible 
at |tii:tp://cxc.harvard.edu/toolkit/pimms.jsp| 
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Figure 2. Optical spectrum of FUTau A from Gemini/GMOS, dominated 
by H« emission and witli several other prominent emission lines indica- 
tive of accretion and/or outflow. The three features between 5300 and 
5500 A are of telluric origin. 



temperature (1.1 keV). The X-ray luminosities of the cool and hot 
components are logLx,i [erg/s] = 29.5 and logLx,2 [erg/s] = 
29.2, respectively, at 0.3 - 8.0 keV. 



2.2 Optical spectroscopy 

FUTau A was observed with Gemini-North and GMOS in long- 
slit mode as part of science program GN-2008A-Q-94, PI Thomas 
Dall. The data were taken in March 2008 and are publicly available 
from the Gemini Science Archive. Two spectra have been obtained, 
one with grism R 600 (exposure time 1200 s, resolution R ~ 4000, 
March 1) and one with grism B 1200 (1200 s, R ~ 4000, March 
16). In both cases, a (/.'5 slit has been used. The quoted values for 
the resolution are taken from the Gemini website. To our knowl- 
edge this source has never been observed wit h a higher spect r al res- 
olution. For comparison, the spectra used bv lLuhman et al.l ( l009al) 
have 7? < 1000 (at 6000 A). 

We downloaded the science and calibration files from the 
Gemini archive and performed standard reduction using the Gem- 
ini/GMOS package in IRAF following the GMOS 'cookbookB 
This includes bias subtraction, flatfield correction, wavelength cal- 
ibration, sky subtraction, and extraction. Since we are only inter- 
ested in emission lines, we did not carry out a flux calibration. 
The final reduced R600 spectrum is shown in Fig. |2] with the most 
prominent emission lines marked. The B 1200 spectrum does not 
contain Ha and has in general lower signal-to-noise ratio. 

The Hq line is broad with an equivalent width (EW) of 146 A, 
which is higher than the values published by Luhma n et al. (009a). 
However, EWs are highly sensitiv e to the chosen valu e for the con- 
tinuum and the integration range dStelzer et aLll2007l) . For a rough 
assessment of the Ha strength in FUTau A and its comparison to 
other objects we calculate the Lna/iboi ratio. We avoid using the 
stellar radius in our estimate because of its uncertainties in the case 
of FU Tau A (see discussion in Sect. |3.5t . Estimating the continuum 
flux from the model and with Fboi = crT'cfft for the bolometric flux, 
we find LHc/Lboi = EW ■ Fcont/Fboi ~ 2 • 10'^ 



http://www.gemini.edu/sciops/data/E^AFdoc/gmosinfospec.html 



We measure a 10 % width jMohantv et alj2005h of 352 km/s 
for the Ha line. The Ha profile is not well-resolved, therefore this 
value has to be seen as upper limit. Still, both the EW and the 10 % 
width are clearly above the typically adopted thresholds between 
non-accretors and accretors of 10 A and 200 km/s. Moreover, FU 
Tau A features a variety of other emission li nes which are often 
seen in accretors, e.g. H/3, He 5876, He 6678 jjavawardhana et al] 
2006). Thus, the optical spectroscopy clearly confirms FUTau A as 
an accreting object. The [01] 6300 and NaD feature may also be 
related to outflow/accretion activity, but contamination by telluric 
emission caimot be excluded. 

From t he 10 % width and based on the correlation given by 
iNatta et alj j2004l) we derived an accretion rate of A^hck — 3.5 ■ 
10~^° M© /yr. Another estimate for the accretion rate can be ob- 
tained ba sed on the line flux in He 5876 and the correlation pub- 
lished by ' Herczeg ~&, Hillenbrand J2008h . In this line, we measure 
EWs of 10 A and 11 A in the two spectra, i.e. no significant vari- 
ability is present. Scaling with the continuum flux, which is ob- 
tained from the A MES-Dusty model s pectrum for T^s ~ 2800 K 
and log 5 = 3.5 jAllard et alj|200lh . this gives a line flux of 
log F5876 [erg/cm'^ /s] = 5.4, which translates into an accre- 
tion rate of Mnei = 7.5 • lO^^^Mo/yr. Adopting Eq. 8 of 
iGullbring et"ai] ( Il998h . the observed A/hcI can be converted into 
the accretion luminosity, iacc.Hci- With an assumed inner disk 
truncation radius i?in = 2 i?, the measured Mho i corresponds to 
iacc,Hci = 1-2 • 10^" erg/s. Our estimate assumes zero veiling in 
the continuum (see Mo hanty et al . 2005); the actual accretion rate 
could therefore be slightly higher. Note, that Mho i is higher than 
Mfia but compatible within the uncertainties given by the empirical 
relations (~ 0.5 logarithmic dex). The agreement with the accre- 
tion rates determined from He 5876 indicates that the large 10 % 
width of the Ha line is intrinsic, i.e. the line is probably resolved. 



3 DISCUSSION 

3.1 X-ray emission from young BDs 

A systematic survey of X-ray emission fro m BDs in the Tauru s 
star-forming region has been performed by iGrosso et al.l ( l2007h . 
They used data from the XMM- Newton Extended Survey of the 
Taurus Molecular Clouds (XEST; iGudel et alj|2007l) that included 
17 young BDs. An X-ray detection fraction of 53 % was ob- 
tained for this sample. Only a small subset of the XEST BDs had 
enough photons collected (> 100 counts) to enable a spectral anal- 
ysis. For most of them the statistics was so poor that the spectral 
shape could be represented adequately with a moderately absorbed 
(iVn < 10^^ cm"^) 1-T thermal model of kT ~ 0.5...1.5 keV, 
typical for the plasma in a magnetically heated corona. 

For none of the BDs detected in XEST the X-ray emission 
is nearly as soft as for FUTau A, despite they all have similar or 
lower absorbing column favoring the detection of a cool plasma 
component. Recall, that in FUTau A the X-ray EM and flux of 
the cooler spectral component exceeds that of the hotter one (see 
Sect. 12.11 ). Dominating soft X-ray plasma has been detected so 
far only in older BDs, e.g. the evolved BD Gl 569 B has an emis- 
sion measure weighted mean temperature of kT — 0.6 keV at an 
age of lOOMyrJ Stelzer 2004) and the 500 Myr-old LP 944-20 had 
0.3 keV kutledge et akboOO) . Both of them were detected during 
a flare, i.e. in an event known to go along with plasma heating, and 
their quiescent emission - if any - is probably even softer. The de- 
crease of I/x/iboi and kT of BDs with age has tentatively been 
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Table 2. X-ray spectral parameters of absorbed 2-T APEC model for FU Tau A compared to those of TW Hya. 



Object 


Instrument 


Xi'ed (dof) 


logA^H 


fcTi 




log EMi 


log EM2 


log Lx* 








[cm^^] 


[keV] 


[keV] 


[cm"-^] 


[cm~3] 


[erg/s] 


FUTauA 


Chandra/ ACIS 


0.9 ( 27) 


21.8-:^ 


n 24" ■ 3"* 
^•^^0.19 


^•^^0.99 


52.9^^;^ 


^'^■•^52.2 


29.5 


TWHyat 


XMM-Newton/EPlC-pn 


2.3 (216) 


20.8 


0.23 


1.22 


53.0 


52.2 


29.8 



* Intrinsic X-ray luminosity in the 0.5 — 2.0 keV passband; distance assumed for TW Hya is 51 pc (Mamajek 2005) 
t The best-fit model of TWHya has peculiar abundances (not listed here); abundances have been set for all elements 
to 0.2 solar for FUTauA because of the low photon statistics. 
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(circles - Giide l et alj i2007|), upward pointin g triangles - IScelsi et"al] 
)2007|)) and BDs (squares- | Grosso et alj JIOOtI) . downward pointing trian- 
gle - iMokler & StelzeJ j2002h ) in Taurus and the position of FU Tau. Non- 
detections in X-rays are indicated by downward pointing arrows. The verti- 
cal lines connect the X-ray emission levels for stars observed and detected 
in two XEST pointings. The luminosities for the TTS given bv lGiidel et alj 
i2007l) and shown in this figure refer to a sHghtly broader energy band 
(0.3 — lOkeV) but for a typical coronal X-ray spectrum (fcT ~ IkeV) 
the difference in flux to the 0.5 — 8.0 keV band used in the BD study is 
negligible. 



explained by reduced coro nal heating related to the increasingly 

cool effective temperatures jStelzer et alj2006l). 

In the 0.5 - 8keV band, used by lOrosso et al.1 ilOOH) for 



XEST, the X-ray luminosity of FUTauA (4 • 10^^ erg/s) is at the 
high end of that of all BDs in Taurus. In Fig. [3] we show the X- 
ray luminosity of the two components in the FU Tau binary com- 
pared to other BDs and TTS in Taurus. Note, that we plot the to- 
tal X-ray luminosity for FUTauA, although below we argue that 
possibly only the weaker, hot component is of coronal origin. If 
only the hotter spectral component of FUTauA, which is more 
similar to the plasma detected from other Taurus BDs, is consid- 
ered FUTauA is still the X-ray brightest BD in Taurus together 
with CFHT-B D-Tau 4. The data for the BDs was extracted from 
1 (l20q7h . those for the TTS are from the catalog of 
( 120070 . In addition, we add data for new Taurus mem- 



Grosso et a 



Giidel et al 



Figure 4. HR diag ram for VLM objects in Taurus including BDs from 
lOrosso e t alj i200 7h and the F U Tau binary on pre-m ain sequence models 
from lBaraffe et alj fT998l) and lChabrier et alj j2000h . X-ray detections are 
shown as filled circles, non-detections are open circles. 



bers identified after the publication of the first results from XEST . 
The X-ray data for these objects are taken from Scelsi et al. ( 20()7|) 
an diMokler & Stelzer ( 200 2). and their stellar parameters are listed 
by iLuhman et alj lo09bl) . Recall, that contrary to some yery-low 
mass (VLM) stars with high X-ray luminosities, indicated by verti- 
cal lines in Fig.|3] FU Tau A did not show variability. 

The fractional X-ray luminosity of FUTauA (logLx/iboi = 
—3.3 in 0.5 — 8keV) is within the range observed for other Tau- 
rus BDs. If only the hot spectral component of FU Tau A is con- 
sidered log (Lx/iboi) = —3.7 which is even more typical for 
yo ung BDs; the medi an for the XEST BDs is —4.0 according 
to iGrosso et alj j2007l) . This 'normal' I/x/iboi ratio despite the 
high Lx is a conseque nce of the par t icular ly high bolometric lu- 
minosity of FU Tau A. iLuhman et al.l ( l009al) give Lboi = O.2L0, 
about a factor 3 higher than for th e optically brightest of the BDs 
examined by Grosso et al, ( l2007h . Its position more than 1 dex 
above the youngest isochrone of the evolutionary calculations by 
lBaraffeetaLHl998l) can be seen in Fig.|4l While several other BDs 
in Taurus are located above the 1 Myr isochrone as well, FUTauA 
is clearly an extreme case. 
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3.2 X-ray emission from substellar twins 

The components in a binary are usually believed to be coeval. For 
BD binaries resolvable with Chandra this enables a study of X-ray 
properties from two substellar objects of different spectral type, i.e. 
mass. In the common hypothesis that the X-ray emission from BDs 
is produced in a corona this means investigating the temperature 
and mass dependence of magnetic activity. 

We detected the optically brighter primary, FUTau A, with 
Chandra at more than two decades higher count rate than the up- 
per limit for the undetected secondary. Considering that these two 
BDs differ by only two spectral subclasses, our result thus implies 
either (i) a very steep decrease of X-ray activity between spectral 
type M7 and M9 or (ii) a crucial influence of some other parameter 
on the efficiency of magnetic activity or (iii) a misinterpretation of 
the link between the two stars, or (iv) an altogether different X-ray 
emission process. Next these hypotheses are investigated. 

Fig.[3] indeed, points at a steep drop in X-ray efficiency around 
Teft ~ 2700 K. However, the decay of goes along with decreas- 
ing optical brightness of BDs, and the sensitivity of XEST and our 
Chandra data for FU Tau do not allow to ascertain whether system- 
atic changes of Lx/iboi occur in the substellar regime at the age 
of Taurus (few Myr). 

Extinction is unlikely the reason for the non-detection of 
FUTauB be cause its Ay is sm aller than that of component A; 
see however IStelzer et al.l |2007) for the case of the BD binary 
2MASS Jl 101 1926-7732383 AB where the component with higher 
Ay was brighter in X-rays than the one with lower A-v- Possi- 
ble parameters that rule magnetic activity are the rotation rate and 
the magnetic field strength and structure. It was shown that pre- 
MS stars are in the saturated regime of th e activity-rotation re- 
lation where Lx i s independent of rotation jPreibisch et"ai]|2005l : 
iBriggs et aLll2007l) . While evolved V LM stars show littl e or no X- 
ray emission despite rapid rotation jBerger etal.ll2010l) . the con- 
nection between X-rays and rotation has not yet been examined 
for young objects in the substellar regime. Rotation periods from 
few hours to several days have been measured for young BDs (e.g. 
IScholz & Eisl6ffelll2004l : |joergens et al.ll2003h . Rotation and mag- 
netic field are unexplored for FUTau; see Sect. l3.5l for an estimate 
of both. 

Another possibility that would explain the strongly different 
X-ray emission levels of FUTau A and B is that the two objects 
are not coeval with the primary being younger and, therefore, more 
active. This interpretation is in line with the HR diagram position 
of FU Tau A above the models while FU Tau B coincides with the 
1 Myr isochrone. An initial age difference of ^5 1 Myr would be- 
come impossible to recognize within a few Myrs. Alternatively, 
FU Tau A and B may be unrelated objects, with component A be- 
ing much closer than the assumed 140 pc. This would explain its 
high X-ray and bolometric luminosity . However, the binary nature 
of this BD system was established bv lLuhman et al.l ( l009ah on the 
basis of the low probability for a chance projection, especially con- 
sidering the high extinction through the B 215 cloud. 

Finally, the different X-ray properties of FUTau A and B 
might point at different emission mechanisms in the two objects. 
In Sect. 13.31 we suggest that a major part of the X-ray emission 
from FU Tau A may be related to the accretion process and not a 
confined corona. 
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Figure 5. Chandra/ ACIS spectrum of FU Tau A and residuals with respect 
to the bestfit presented in Table |2] (red in color print) compared to the 2-T 
bestfit model for TWHya (solid line; blue in color print). The spectrum of 
TWHya has been scaled by a factor of 0.07 to make the comparison with 
FU Tau A easier. We also have artificially increased the absorbing column 
density of TW Hya to the value measured for FU Tau A and reduced the 
neon abundance from its atypically high value to solar. 



3.3 FU Tau A, a substellar analog to TW Hya? 

In Fig.|5]we display the X-ray spectrum of FU Tau A. As mentioned 
already in Sect. 2.1, the spectral fitting revealed an unexpectedly 
strong soft plasma component. This characteristic is reminiscent of 
the TTS TW Hya. In Table |2] we compare the best-fitting param- 
eters of FU Tau A to those obtained from an XMM-Newton/EPlC 
observation for TWHya. Despite the best fit result presented by 
[Stelzer & Schmitt (2004) for TWHya was a 3-T model, we show 
here the parameters for a 2-T model for better comparison with 
FU Tau A. This model spectrum for TW Hya is overlaid in Fig. |5] 
on the observed spectrum of FU Tau A. (A direct graphical compar- 
ison of the two observations is impossible because of the different 
instruments with which the data were obtained.) To ease the com- 
parison, the model of TWHya was multiplied with a factor 0.07 
(accounting for the difference in distance and the factor two dif- 
ference in X-ray luminosity), the column density was increased to 
the value observed for FU Tau A, and the neon abundance was de- 
creased to the solar level. There is almost a 1 : 1 correspondence 
between the spectrum of FUTau A and the down-scaled, artificially 
absorbed spectrum of TW Hya, demonstrating that the temperature 
structure is very similar. 

Both objects are dominated by soft emission with only a weak 
emission component from the higher energy photons that are typ- 
ical for stellar coronae. By means of high-resolution X-ray spec- 
troscopy, plasma densities (Uc ~ 10^^ cm^^) two dex hig her than 
in ste llar coronae have b een m easured for TW Hya 1 Kastn er et al.l 
I2OO2I ; IStelzer&Schniittll2004h . Its soft X-ray emission and the 
measured high densities are compatible with the conditions in an 
accretion shock. This way TW Hya was established as the prototype 
of a pre-MS star with X-ray emission from accretion columns. A 
handful of other accreting TTSs have since be en recognized to dis- 
play high density X -ray emitting plasmas (e.g. lSchmitt et al.ll2005l ; 
iGunther et ai]|2006l) , and an excess of soft X-rays over the typi- 
cal coronal emission of young stars has been detected in a num- 
ber of accreting TTSs l lGUdel & Telleschl2007l) . Obtaining a high- 
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resolution X-ray spectrum for FUTau A or any other BD is out-of- 
reach for present-day X-ray instrumentation, such that its plasma 
density can not be constrained. However, on the basis of the similar 
temperature structure between the X-ray spectrum of FU Tau A and 
TW Hya we speculate that the soft X-ray emission from FU Tau A, 
represented by the dominating component # 1 in the spectral fit, 
may be produced in accretion shocks rather than in a corona. 

To check if the accretion scenario is compatible with the ob- 
served properties we compute the theoretically expected emission 
measure from the post-shock region, 

£Mp,h = Afacc ■ ^^^p,^"^'' , (1) 

^mpi-'(ipsh) 

where fca is the Boltzmann constant, fj, — 0.61 the mean parti- 
cle weight in units of the proton mass m p, and P(Tpsh) th e ra- 
diative loss function extracted from APEC jSmith et alj200lh . The 
post-shock temperature is Tp^h ~ ^ ^i^'^O' where vo is the infall 
speed. For the accretion rate. Mace, we substitute the value ob- 
tained from He 5876, MhcI. Adopting for no the velocities mea- 
sured from the width of the Ha profile, vua ~ 175km/s, we find 
fcTpsh ~ 0.03 keV and log EMpsh [cm"^] = 52.7. While the pre- 
dicted emission measure is roughly comparable to the soft X-ray 
plasma, EMi from Table (2] the post-shock temperature is smaller 
than the observed X-ray temperature, kTi. The thermal emission 
connected to the low temperature expected from material infalling 
towards a BD should, indeed, be too soft to be significant in the 
X-ray band. On the other hand, as discussed in Sect. l3. li the domi- 
nance of the soft X-ray component is completely untypical for mag- 
netic activity on such ayoung BD. So, if accretion is not responsible 
for FU Tau A's X-ray emission, it must have a very unusual corona. 



3.4 Accretion properties of FU Tau A 

FUTau A stands out among BDs of similar mass for its high 
accretion rate. The empirical M — Afacc relation has a large 
spread of controversial origin. From optical spectroscopy, values 
of -Mace ~ ]^Q-i2...-9 y[^/yj- have bee n measured for BDs with 
M ~ 0.05 M© (e.g. lNattaetZll2006h . However, in that mass 
range accretors with Mace ^ lO~^''M0/yr are fo und only in 
pOph , the youngest of the examined star-forming sites jNatta et al] 
l2004h . and the M* - Maec relation predicts log Aface [Mo/yr] = 
— 10.7, more than 1 dex below our value for A/hc. 

Strong accretion produces hot spots on the surface of the star 
that might introduce brightness variations in the course of a rotation 
cycle. The two epochs of Sloan Digital Sky S urvey (SDSS) photom - 
etry (from Dec 6 and 29 2002) presented bv lLuhman etai] j009ah . 
indeed, indicate strong flux variations. The amplitude is 0.07 mag 
in the z band and increases monotonously towards blue wave- 
lengths (e.g., 0.27 mag in r band, 0.44 mag in g ba nd). The wave- 
length dependence from a standard extinction law jMathisiri990l) 
is not steep enough to explain this behaviour with variable extinc- 
tion and we assume that we see the signature of star spot s. U sing a 
simp le spot model comparable to the one presented bv lScholz et al.l 
( 12005 ) and assuming blackbody spectra for the spotted and unspot- 
ted surface area, we find that cool spots with AT = 0— 1000 K and 
filling factor / = 5 — 50 % do not provide the observed steep de- 
crease of amplitude with wavelength. On the other hand, hot spots 
(caused by the accretion flow) with AT > 500 K and / = 5 - 10 % 
yield a decent match to the observed spectral dependence of the 
amplitudes. We assume here that the minimum photometry corre- 
sponds to the photosphere which may not be the case. Clearly, the 



two epochs of SDSS data provide only a lower limit to the vari- 
ability of FUTau A. Spot s with AT of a few 1000 K, as have been 
measured on other BDs / Scholz et alj|2009 '). can not be ruled out 
for FU Tau A on the basis of existing photometry. Hot spots with 
these characteristics have a luminosity of Lspot = 1-5 ■ 10'^'^ erg/s, 
i.e ~ 20 % of the bolometric luminosity of FU Tau A. 



3.5 The relation between stellar parameters and 
accretion/activity in FU Tau A 

One puzzling feature of FUTau A is its position high above the 
youngest isochrone in the HR diagram. An accurate assessment 
of (sub)stellar parameters is crucial for the evaluation of accretion 
and X-ray properties as can be seen from the following open prob- 
lems. First, the velocities derived in Sect. 13.31 from the Ha line 
are difficult to reconcile with the free-fall speed for the published 
(sub)stellar parameters of FU Tau A. In particular its large Lboi and 
ensuing huge radius (1.8 Rq) imply vg ~ lOOkm/s, almost a fac- 
tor two lower than the value derived from Ha. Secondly, the spot 
luminosity and the luminosity derived from the kinetic energy of 
the infalling material (Lace.Hc i) could be better reconciled assum- 
ing a 1 dex fainter Lboi for FUTau A; this would moreover make 
the BD compatible with the 1 Myr isochrone. On the other hand, 
the high bolometric luminosity of FU Tau A comes along with high 
X- ray luminosity, supportin g an unexplained trend first noticed 
by IStelzer & Micelal ( l2007h who found that most of the BDs in 
star-forming regions detected in X-rays are located on the 1 Myr 
isochrone or above the evolutionary calculations. This trend is also 
seen in Fig.|4] Thus, FUTau A may not be an isolated albeit an ex- 
treme case. In the following we consider various possible causes 
for the high luminosity of FU Tau A: magnetic activity, accretion, 
problems with models and wrong distance. 



3.5.1 Reduced convection affecting temperature and radius 

A mechanism that can produce larger radius (and cooler tempera- 
ture) for a given mass, without affecting the bolometric luminosity, 
i.e. shifting the object horizontally in the HR diagram, is reduced 
convection efficiency and ensuing heat flux because of fast rotation 
an d/or strong magnetic fi elds. Su ch a scenario has been in vestigated 
bv lChabrier et all ( l2007h and bv lMacDonald & MullaiJ {2009) us- 
ing different approaches to model the influence of the magnetic 
field. Both studies have successfully explained the temperature re- 
versal of the eclipsin g spectroscopic binary 2M 05352184-0546085 
dStassun et al]|2007h with strong magnetic effects on the primary. 
Indeed, in that object the primary which has a lower effective tem- 
perature than the secondary is strongly active in terms of h igh rota- 
tion velocity and strong Ha emission l iReiners et al.ll2007l) . 

Shifting FUTau A to higher temperatures onto the IMyr 
isochrone in the HR diagram to make up for a presumed effect of 
activity on T^s and i?, yields a 'true' mass of ~ 0.2 Mq and 'nom- 
inal' temperature of ~ 3190 K. We can estimate the rotation period 
and field strength required for magnetospheric accretion models. 
We adopt the X-wind model that defines the disk truncation ra- 
dius J?x as the distance from the star where the Keplerian angu- 
lar velocity equals the stellar angular velocity (see lOstriker & Shti 
ri995; Mohanty &Shu 2008). Assuming i?in = 2 R, we obtain 
with the radius of 1.83 -R© derived from luminosity and temper- 
ature and with the 'corrected' stellar mass of FUTau A a period 
of Prot ~ 1.8 d. Using the model of lKoenigll h99lh where the 
disk truncation radius is smaller than the corrotation radius yields 
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a longer period, Prot ~ 5 d. We note in passing that a larger in- 
ner disk radius (i?in = 5 /?, is often used in the literature) yields 
periods between 7. ..21 d depending on the details of the accretion 
model. These latter values for the period are at the high end of the 



r^l 



values measured for BDs and TTS (e.g . Scholz & Eisl6ffelll2004l : 
iHerbst et al.l2007l) . Indeed. iRebull et alj ^20061) presented evidence 
that stars with disks have longer periods on a statistical basis than 
stars without disks but this does not exclude the exi stence of (s ome) 
fast-rotating accretors. For the observed Mhc i . iKoenigll ( 1 199 ll) pre- 
dicts for FUTauA a magnetic field between B ~ 85... 400 G, 
for assumptions of iiin = 2 R, and 5Rq, respectively. Mag- 
netic fields with kG strength are not unusual for TT S ( Johns-Kri^ 
12007 ). For BDs with ages of ~ 1 - lOMyr Reiners et akj <200' 
could not find detectable fields with upper limits of ^ 1 kG for all 
of the accretors, while young but non-accreting BDs and evolved 
ultracool dwarfs in the field exhibit field strengths of a few kG 
( iReiners & Basri 2007). Howe ver, the sample of only three accre- 
tors studied by IReiners et alj |2009[) is too small to be clear evi- 
dence against the presence of ^5 1 kG fields on accreting substellar 
objects. 



To summarize, for a reasonable range of _Rin we obtain ei- 
ther fast rotation or strong magnetic field for FUTauA, i.e. at least 
one of the two conditions for suppressing convection seems to be 
fullfiUed independent of the position of the disk truncation radius. 
The high level of (X-ray) activity of FU Tau A - even neglecting 
the dominating cool X-ray plasma possibly produced by accretion 
- is compatible with this scenario. Similarly, the high value for the 
Ha flux corresponds to a high Lua/Lhoi ratio. However, given the 
dominant contribution of accretion in forming the Ha line this mea- 
surement can not be used as a diagnostic for strong magnetic activ- 
ity in FU Tau A, contrary to the case of the non-accreting eclipsing 
binary 2MASS 05352184-0546085 (Reiners et al. 2007). Cool star 
spots are expected as a manifestation of inhibited convection. As ar- 
gued in Sect. 13.41 the two epochs of photometric observations can 
not be explained by cool spots alone, but they can not rule out the 
presence of cool spots either. More detailed monitoring is required 
to test this scenario. 



While the above considerations are not conclusive, in the ab- 
sence of measured rotation, field strength and good constraints on 
the nature of the spots it is not excluded that rotation/fields affect 
the stellar parameters of FU Tau A. However, the presence of accre- 
tion, contrary to the case of 2M 05352184-0546085, certainly com- 
plicates the situation. A quantitative theoretical evaluation of mag- 
netic field effects in this particular object would certainly provide 
further insight. We now discuss a few implications of this scenario 
for the interpretation of our X-ray and Ha observations. First, the 
higher mass would move FU Tau A into the stellar regime. Higher 
mass objects have higher X-ray luminosity, i.e. FUTauA would 
be less extreme in terms of X-ray brightness. The I/x/iboi ratio 
would be unchanged and it is, indeed, compatible with the values 
observed in young VLM stars; in the Orion Nebula Cluster the me- 
dian of 0.2 Mq stars is log ( Lx/Lboi) ~ —3.5 an d the spread is 
at least one logarithmic dex jPreibisch et alj[2005l) . Secondly, the 
free-fall velocity would be higher (~ 200 km/s) and in better agree- 
ment with the observed Ha profile. Moreover, the higher velocity 
corresponds to a predicted post-shock temperature which is more 
similar to the X-ray observed temperature. Finally, the accretion 
luminosity derived from Mhb would increase by a factor four but 
still be much smaller than the estimated spot luminosity. 



3.5.2 Excess accretion luminosity 

Alternatively, the position of in the HR diagram for FU Tau A and 
simi lar objects may be du e to excess luminosity related to accre- 
tion. iLuhmanetalltoah have estimated Lboi for FU Tau A from 
its J band magnitude, where deviations from photospheric emis- 
sion due to the presence of a disk and accretion are usually small- 
est. However, strong (accretion) variability in the J band can not 
be ruled out with the available data. 

If the true bolometric luminosity of FU Tau A is significantly 
lower than the adopted value the radius is smaller (for given mass 
of 0.05 M0) and the free-fall velocity gets larger, e.g. a factor 
10 smaller Lboi gives vg = 180 km/s in perfect agreement with 
the observed Ha width. On the other hand, this interpretation 
implies a high fractional X-ray luminosity for FUTauA and the 
other BDs that are located above the models in the HR diagram 
(ix/iboi ^ 10"^). Such high X -ray emission levels are difficult 
to explain in terms of T Tauri like activity. Moreover, two of the 
three Taurus BDs that are closest in the HR diagram to FU Tau A, 
i.e. highest above the models, are not accreting based on the equiv- 
alent width of their Ha emission ( Grosso et al.l2007h and they even 



have no disks according to lLuhman et al.l ( |2010|) . making accretion 
luminosity an unlikely explanation for their discrepancy with the 
evolutionary models. 

3.5.3 Shortcomings of evolutionary models 

It is well-known that different sets of pre-main sequence models 
are inconsistent among each other and none of them describes ade- 
quately all observations. It is not clear to what extent these discrep- 
ancies between data and model are due to problems with the stellar 
param eters derived from ob servations or due to shortcomings in the 
theorv. lBaraffe et al] j2009l) showed that the typically observed age 
spread of 1 — 10 Myr in the HR diagram for stellar populations in 
star forming regions can be explained by episodic accretion dur- 
ing the protostellar phase. However, they can not explain the pres- 
ence of objects above the 1 Myr isochrone. These authors suggest 
that such objects are very young, having experienced their episodic 
protostellar accretion events quite recently and not yet contracted 
to the 1 Myr position of non-accreting objects of the same mass. In 
this view the two components of the FU Tau binary are not coeval. 

An age difference of a few hundred thousand years has 
been detecte d for the first t ime in an equal-mass eclipsing binary. 
Par 1802, by IStassun et al.l ( |2008|) . demonstrating that subsequent 
formation of the components in stellar systems is possible. For the 
case of FU Tau the existing evolutionary models are not helpful for 
constraining the possible age difference. Compared to Par 1802, 
FUTau is of lower mass and much wider separation. Moreover, 
the masses of the two components are not equal introducing addi- 
tional uncertainties in an evalution of their evolution. To our knowl- 
edge, star formation theories have not yet made predictions on such 
young non-coeval binaries, their properties and possible influence 
of the environmental conditions. 



3.5.4 Wrong distance or binarity ofFU Tau A 

Finally, we mention for completeness further scenarios that would 
reduce the bolometric luminosity of FUTauA and bring it in ac- 
cordance with the 1 Myr isochrone: (i) A much closer distance than 
assumed; this would also result in a lower, and more typical. X-ray 
luminosity. However, the ensuing non-binarity with component B 
is countered by statistical arguments, and distance is unlikely to be 
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responsible for similar cases of other BDs above the evolutionary 
models, (ii) The primary in the FUTau system might be a binary; 
however, this would not provide a solution because it would de- 
crease Lboi by at most a factor of two. 



4 SUMMARY 

Among the most puzzling features in the enigmatic BD FU Tau A 
are a very high bolometric luminosity for its effective temperature 
and its strong and unexpectedly soft X-ray luminosity pointing at 
an origin in accretion rather than a corona. 

Assuming a mis-interpretation of the stellar parameters as dis- 
cussed in Sect. 13.51 would resolve many open problems: HR dia- 
gram position, free-fall versus observed velocities, spectroscopic 
and photometric accretion luminosities. Suppressed convection as 
a result of fast rotation and/or strong magnetic field might be re- 
sponsible for the apparent excess luminosity of FU Tau A with re- 
spect to pre-main sequence evolution calculations. If confirmed this 
would be the first case of an accreting BD with strong field effects. 
However, in the absence of observational constraints on the rota- 
tion rate and field strength this explanation is not much more than 
a plausible hypothesis. Keeping in mind the many pieces of evi- 
dence for strong accretion in FUTau A from optical photometry, 
optical spectroscopy and X-ray emission, it is possible that the ex- 
cess luminosity may be explained by a significant contribution of 
accretion to the J magnitude. Further constraints on the accretion 
properties of this benchmark BD are required to give our estimates 
a quantitatively sound basis. Last but not least, the existence of 
other BDs that have stellar parameters incompatible with evolution- 
ary models, albeit less pronounced, may suggest problems with the 
ages predicted by evolutionary models. If FU Tau A is younger than 
1 Myr and not coeval with FUTauB, as suggested by Baraffe et al.l 
l l2009l) , its extraordinarily strong and soft X-ray emission, irrespec- 
tive of whether due to accretion or to activity, could then be a result 
of its young age. In the near future we aim at further constraining 
observationally the accretion and activity of this unusual object. 

Our Chandra observations of FUTau A represent the first in- 
dications for X-ray emission from accretion shocks in a BD (pro- 
vided FUTau A is a substellar object and not a star as argued 
in Sect. 13.5.1b . We note, that even if free-fall velocity and Ho? 
width of FU Tau A can be reconciled with the scenarios described 
in Sect. 13.51 the predicted X-ray post-shock temperature is still 
lower than our observed value. While this discrepancy remains un- 
explained, a similar case is represented by the TTS Hen 3-600. 
This object shows clear evidence of accretion from its Ha char- 
acteristics as well as the density and s oftness of its X-ray emit- 
ting plasma dHuenemoerder et"^l2007h . However, quantitatively 
its mass (0.2 Mq) and radius {R = 0.9 Rq) imply a post-shock 
temperature of 1.1 MK, to be compared to the measured X-ray 
plasma temperature of ~ 3 MK. 
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